Search for: All records

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering — removing sequencing bases, reads, genetic variants and/or individuals from a dataset — to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy–Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima’s D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne). 

Abstract Divergence is often ephemeral, and populations that diverge in response to regional topographic and climatic factors may not remain reproductively isolated when they come into secondary contact. We investigated the geographical structure and evolutionary history of population divergence withinSceloporus occidentalis(western fence lizard), a habitat generalist with a broad distribution that spans the major biogeographical regions of Western North America. We used double digest RAD sequencing to infer population structure, phylogeny and demography. Population genetic structure is hierarchical and geographically structured with evidence for gene flow between biogeographical regions. Consistent with the isolation–expansion model of divergence during Quaternary glacial–interglacial cycles, gene flow and secondary contact are supported as important processes explaining the demographic histories of populations. Although populations may have diverged as they spread northward in a ring‐like manner around the Sierra Nevada and southern Cascade Ranges, there is strong evidence for gene flow among populations at the northern terminus of the ring. We propose the concept of an “ephemeral ring species” and contrastS. occidentaliswith the classic North American ring species,Ensatina eschscholtzii. Contrary to expectations of lower genetic diversity at northern latitudes following post‐Quaternary‐glaciation expansion, the ephemeral nature of divergence inS. occidentalishas produced centres of high genetic diversity for different reasons in the south (long‐term stability) vs. the north (secondary contact).